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SUMMARY 


The flow of fluids with and without heat sources and subject to 
body forces between two plane parallel surfaces which are oriented in 
the direction of the generating body force is analyzed under the condi- 
tion that the temperature vary linearly along these surfaces. It is 
found that a modified Rayleigh number (product of the reciprocal of the 
ratio of specific heats and the Prandtl number Pr and the mo di fied 

. pfyd 

Grashof number Gr/J an well as a parameter = Pr Gr^ is of 

C P 

significance in this problem; where p is the volumetric expansion 
coefficient, fx is the negative of the X-component of body force per 
unit mass, d is the characteristic length, and Cp is the specific 
heat at constant pressure. Solutions of this problem are obtained in 
terms of "universal" functions which are tabulated for simple applica- 
tion to specific cases. Representative velocity and temperature dis- 
tributions from which detailed study of the heat transfer is made are 
then computed. When the ratio of CK^ (where C is related to the mass 
flow) to the Rayleigh number is of unit order of magnitude, the effects 
of aerodynamic or frictional heating can be appreciable. Asymptotic 
solutions (for large values of the Rayleigh number) which render the 
computations simple are also presented. 


Comparison of the results from the method given herein with those 
obtained elsewhere in an approximate manner for a special case simu- 
lating the natural -convection flow of fluids with heat sources in a com- 
pletely enclosed region shows that the approximate method is suffi- 
ciently accurate for problems in which the modified Rayleigh number is 
less than 10^. 
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INTRODUCTION 

In recent years the transfer of heat to and from enclosed or par- * 

tially enclosed regions "by means of natural convection or "by a combina- 
tion of natural and forced convection has taken on new significance in 
the fields of aeronautics, atomic power, electronics, and chemical en- 
gineering. Most of the information on these modes of heat transfer under 
such conditions is of a semiempirical or specialized nature; relatively 
little detailed information exists for internal natural-convection flows. 

In reference 1 there appears one of the few attempts to determine theo- 
retically the velocity and temperature distributions in detail and hence 
the heat transfer for an internal flow problem of this kind. In that 
reference a solution was found for the fully developed flow of fluids 
with and without heat sources between two long parallel plates with con- 
stant wall temperatures (where one could be different from the other) 
oriented in the direction of the generating body force. The information 
obtained therein is of practical value in connection with fully developed 
flows subject to body forces where the surfaces are maintained at uniform 
temperatures. Lighthill (ref. 2) employs integral methods to study the 
natural -convection flow in tubes with either one end or both ends closed 
and with constant wall temperatures, and in reference 3 approximate 
superimposed free- and forced-convection flows are obtained for short 
channels and pipes. 4 

As the next step in the study of natural convection or combined 
natural and forced flows in confined spaces, consideration Is here given 1 

to the configuration of reference 1 with the exception that the thermal 
boundary condition specified is that the surface temperatures vary 
linearly along the plates or surfaces. (One surface, however, may be at 
a different local temperature from the other but the slopes of the tem- 
perature distributions on each surface are taken to be equal.) The 
analogous forced-convection problem is treated in references 4, 5, and 6. 

The present problem simulates several important physical occurrences of 
this phenomenon; for example, it could represent the case where the out- 
side of the channel formed by the plates Is cooled (or heated) by a 
counterflow. In addition, the present problem represents a more general 
case than was considered in reference 1, since here the temperature will 
no longei be restricted to be a function of the transverse coordinate alone 
and, hence, energy convective as well as mass convective effects will be 
included. 

The solution is obtained in terms of functions which depend on only 
one of the several associated dimensionless parameters, and these func- 
tions are tabulated so that specific cases can be easily computed. Solu- 
tions for pure natural convection and for superimposed natural and forced 
convection are shown to be essentially Identical. Representative velocity 
and temperature distributions are also presented, and the effects of 
frictional or aerodynamic heating on the flow and heat transfer are 
discussed. 
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A special case which simulates a completely enclosed region in which 
there is no net mass flow, the vails are at the same temperature, and 
heat is generated uni formly "by heat sources is treated, and detailed 
velocity and temperature profiles are obtained from which the heat trans- 
fer is determined. This special case was treated in an approximate 
manner in reference 7. 

Consideration is also given to the problem of convective inversion, 
that is, to the cases where the modified Grashof number changes sign. 

It is shown that convective inversion due to changes in the body force 
direction, to changes in the sign of the volumetric expansion coeffi- 
cient, or to changes in the sign of the axial temperature gradient alters 
the character of the problem, because, it is believed, under these con- 
ditions the flow becomes unstable because of heating from below (see 
pp. 104 to 107, ref. 8). 


ANALYSIS 

Formulation of the Problem 

* The study to be made here is that of the laminar fully developed 

flow of fluids with and without heat sources and subject to a body force 
between two plane parallel surfaces open at both ends and oriented in 
the direction of the generating body force (see fig. 1). It Is further 
specified that there shall be linear (with equal slopes) temperature 
variations along the walls but that the walls need not necessarily be at 
the same temperature. The flow is assumed to be parallel to the axis of 
the channel (that Is, the only nonvanishing velocity component is the one 
in the longitudinal direction) and in addition it Is assumed that the 
physical properties (for example, c p and p) of the fluids are constants 
and that the essential influence of the density changes on the flow is 
taken into account by the introduction of the volumetric expansion co- 
efficient in the body force term (that Is, the other influences of vari- 
able density and the variation of the expansion coefficient with tempera- 
ture are negligible) . Discussions of the justification of the assump- 
tions can be found in references 1, 9, and 10. 

Under the conditions stated, the basic equations with body forces 
included expressing the conservation of mass, momentum, and energy (see 
ref. l) become, respectively, 


ax 2 



\ 


( 2 ) 
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5P 

Sy 


0 


(3) 


d 2 T* d 2 !* P°v Tr dT* u /2u\ 2 Q 
dX 2 dY 2 * k U " k \5Y J " k 


(4) 


(See appendix A for a complete list of the symbols used herein. ) With 
cognizance taken of equations (l) and (3) the above system reduces to 


d% 

dY 2 



(5) 


5%* * 5 2 t* p°v St * u /du\ 2 

SX 2 SY 2 “ k ST - & \dYy 

where the velocity is a function of the transverse 
and the pressure is a function of only the longitudinal coordinate. 

The body force term in equation (5) can be written as a buoyancy 
term by introducing the volumetric expansion coefficient p in a manner 
similar to that described in appendix B of reference 1. Equation (5) 
then becomes 


Q 

k 


( 6 ) 


coordinate Y only 


* 


§ ♦ ^ - Tw 0 ‘) “!(I + (7) 

where the subscript w 0 refers to the surface at Y = 0 (see fig. l) . 

The boundary conditions associated with this problem are that the 
velocity at the walls must vanish (the no-slip condition for viscous 
fluids) and that the temperature must vary linearly along the walls. 
(Note that the latter condition implies that the temperature gradients 


along the walls and hence the axial heat flux along the walls must be 
constant.) In order to satisfy the temperature conditions and equation 
(7 ), the temperature must be of the form 


T*(X,Y) - AX + T(Y) 

(8) 

Mathematically, 

the boundary conditions are formulated as 



U(0) a u(d) = 0 

(9) 


T*(x,0) = AX + T(0) 

(10) 


> 
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T*(X,d) =» AX + T(d) (11) 

Substituting equation (8) into equations (7) and (6), respectively, 
yields 


and 


where 


d% + PP f X 

dY 2 + ^ 


5 (S + V*) 


( 12 ) 



PCyA 


u + 




k VdX 



(13) 


0 * T(Y) - T(0) 


Since the left side of equation (12) is a function of Y alone 
and the right side is a function of X alone, it is clear that each 
side must be equal to a constant. Thus, equation (12) can be written as 


d % 
dY 2 



or 


(14a) 


J (§ + P» 0 JS) - 0 (I4,b) 

The temperature now appears in equations (13) and (14) as 6 which 
is independent of the longitudinal coordinate X, and hence the thermal 
boundary conditions can be written as 


0 ( 0 ) = 0 

0(d) = T(d) - T(0) = 0 Wl (15) 


To nondimens ionalize equations (14) and (13), 


let 
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U 


0 


ku l c p AA 

pPf x d 2 = U ^A 

Y - yd 

_ Ad 

p 2p2 fx 2 d 4 K a 




> 


J 


(16) 


where K A 


PrGr A 



is the new dimensionless parameter presented 


in reference 1, hut here the Grashof number is modified because it is 
based on Ad rather than on a temperature difference , that is, 

Gr A - Hence, equations (14), (13), (9), and (15) become 


u” + T = CK A (17) 

T" -Eau+ (u') 2 + oKa = 0 ( 1S ) 

u(0) =* u(l) = 0 (19) 


T(0) « 0 


( 20 ) 


7(1) 


K A Svi 

Ad 


( 21 ) 


where the primes denote differentiation with respect to y, 

Ra = (l/r) FrGr A is -a modified Rayleigh criterion or number (see p. 105, 

re f. 8), a, a Qd/kA is the dimensionless heat source parameter, and 

C = C /^Prd 3 /cpAK A . Eliminating T between equations (17 ) and (18 ) 

results in 


u lv - (u 1 )^ + Ra u - aK A b» 0 


( 22 ) 
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with the boundary conditions 


u(0) = u(l) = 0 

(23) 

u" (0) = CK a 

(24) 

u"(l) = (C - ®wi/Ad) K A = mCK A 

(25) 


where m = 1 - 0 w /CAd. The constant C which appears as a parameter in 
the "boundary-value problem described by equations (22) to (25) merely 
specifies the temperature level (see eq. (17)) of the problem. In order 
to define completely the temperatures and velocities, this constant must 
be related in some way to the physics of the problem. From equation (l4b) 
it can be seen that C could be determined from the pressure gradient 
along the channel; that is, C is essentially connected with the end 
conditions to which the channel is subject. Since the pressure gradient 
may not be known a priori, in the subsequent section dealing with the 
solution of the present problem C will be related to the end conditions 
by the mass flow in the channel, which re main s invarient over the entire 
length of the channel. 

Note that solution of the preceding boundary-value problem will 
yield velocity and temperature distributions for both natural- convection 
and combined natural- and forced-convection flows . The forced- 
convection pressure gradient merely alters the magnitude of the constant 
C. A discussion of such a superimposed flow problem under special con- 
ditions is given in reference 11. 

Several interesting observations can now be made concerning these 
equations. First, comparison of equation (22) with the corresponding 
equation in reference 1 shows that they are Identical except for the 
third term in equation (22), which does not appear at all In the equa- 
tion in reference 1. This term stems from the convection term in the 
energy equation; hence, in the present problem energy convection effects 
will be Included. The energy convection term vanished identically in 
reference 1 because of the assumption that the velocity and temperature 
profiles were independent of the axial coordinate. Second, since the 
convection term appears with a coefficient, another dimensionless param- 
eter (Ra, the modified Rayleigh number) is associated with this problem 
and its influence on the results must be studied. Finally, the condi- 
tions of the problem require that the temperature be of the form 
T* = AX + T(y) (see eq. (8)), and hence the longitudinal heat flux Is 
everywhere constant. 
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The bo undar y- value problem stated In equations (22) to (25) can be 
written in more convenient form by defining 

v =» Ra u - a K A (26) 

» 

Hence 

-iv _ 1 , 2 , . 

v + Ra v - (v<) =0 (27) 

v(0) = 7(1) = - ccK a (28) 

7"(0) = + Ra CK a v"(l) = + m Ra CK a (29) 


Solutions of the Boundary-Value Problem 

Equation (27) is nonlinear (the nonlinear term is due to the fric- 
tional or aerodynamic heating) and therefore, as in reference 1, a method 
of successive approximations will be employed to find its solution. To 
this end, equation (27) is written 


_iv 


+ Ra C^A-i) 2 - 0 


(30) 


where n => 0,- 1 denotes the particular term in the approximation 


Vq + V]_ and v'^ = 0. Let 


1 - 
64 V 


2y - 1 


oK fl 


X = 


J = 


64 

Ra 1 / ^CKi 


v _ Ra 
R “ 16 




( 31 ) 


J 


64 
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Equations (30), (28), and (29) "become, respectively, 

v£ V + Ev n - i(v i . 1 )2 = 0 (32) 

V 0 (-1) = V 0 (l) = X ( 33 ) 

vg(-l) - J \fR Vq ( i) = mJ*jR (34a) 

V n (-1) * V n (l) = v"(-l) = v”(l) = 0 n^O (34b) 


where the subscripts now denote differentiation with respect to tj. 

Zeroth-order approximation. - In the zeroth-order approximation, 
the nonlinear term which is associated with the frictional heating does 
not appear and the problem then consists of solving the equation 

vj v + Ev 0 = 0 (35) 

subject to the boundary conditions 

v Q (-l) = v 0 (l) = A (36) 

vg(-l) = vj(l) = mJ^R (37) 

For simplicity of computations the solution can be obtained in terms of 
symmetric and antisymmetric functions of b, depending only on the param- 
eter R (or Ra), by setting 

(38) 

Vq j (where 



where the boundary conditions to be satisfied by the various 
j = 0, 1, 2) are given in the following table: 


Subscript v(-l) 

v(l) 

v" (-1) 

V" (1) 

00 1 

1 

0 

0 

01 0 

0 

-#/F 

a/r 

02 0 

0 


a/F 
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The functions v^ must each satisfy equation (35), which has the general 
solution for positive values of the Rayleigh number 

Vq = a x exp(errj) + a 2 exp(-erTj) + a 3 exp(grt]) + a 4 exp(-gni) (39) 

where 


e =* i + 1 
S - i - 1 


and 

r = R 1 / 4 /^ 

Using the boundary conditions to evaluate the constants in equation (39) 
and expressing the solution in terms of real products of circular and 
hyperbolic functions show that 

a 5 i = (cosh r cos r cosh rT] cob rT] 

cosh 2 r + cos 2 r - 1 

+ sinh r sin r sinh ri} sin rq) (40) 


cos^r - cosher 


2 ~ (cosh r sin r sinh rq cos rq 


si nh r cos r cosh rT} sin rr]) 


(41) 


v nC) = o o (cosh r cos r sinh ri] sin rri 

cosh r + cos r - 1 

- sinh r Bin r cosh rT] cos rr]) (42) 

Note that Vqq and Vq 2 are symmetric functions and Vq^ is an anti- 
symmetric function of T], These "universal" functions v Q j are given 
for various values of Ra in table I. 
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First-order approximation . - In order to include the effects of 
frictional heating to a first order , let n = 1 in equation (32). Thus 


v* v + Rv 1 = | (v «) 2 = f ( rj ) (43) 

where, in view of equation (38), f (t}) is given "by 

f (il-) = | [X V + XJ(m + l)y^ v* 2 + v-^2 

j2(m + 1)2 2 + _ 3W1 v i + v-t v» 

4 02 k i;v 00 01 2 01 02 

(44) 


The "boundary conditions on v^ are 

V-l(-I) = v x (l) - v"(-l) = v£(l) = 0 
Once again, to obtain "universal-type" functions let 

v 1 = A 2 v 10 + X J(m + l)^ + 2 v i2 


(45) 


. J 2 (m + l ) 2 . , w , x . J 2 (m 2 - 1) 

+ . - S - - J Vl3 + *j(m - l)v M + 5 i 


'15 


(46) 


In view of the form of equations (43) to (46), the functions v 10 , v-, , 
v-^g, and v^g are symmetric; and v ^ 4 and v^g are antisymmetric 
functions. Each of the functions v-^ (k = 0, 1, ... 5) will "be a 

solution of equation (43), "but where the nonhomogeneous term is only 
the related part of equation (44) (for example, the first term on the 
right-hand side of equation (44) is associated with the v^q solution 

and so forth) . 

A particular solution of equation (43) can "be constructed from 


( V lk)p f k(£> G ( T 1 - £ ) d £ 
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where G is the Green's function which satisfies the homogeneous part of 
equation (43) and the conditions G(0) « G* (0) =* G" (0) «■ 0, G* * * (0) « 1. 
These arbitrary conditions were applied to yield a simple form for the 
Green's function. Hence 


(v ) 

*■ lkP 



sinh g r(t} - I) - g sinh e r(T] - |)J f fc (£ )d| 


(47) 


Note that if the f fc is a symmetric function, (v^)p is also symmetric; 

and if fjj is antisymmetric, so is (yLk)p* From the boundary conditions, 
equation (45), it can be seen that v i > must be either symmetric or anti- 
symmetric. Therefore, if f^ is symmetric, the complementary solution 
used with the particular solution must be symmetric; and if f^ is anti- 
symmetric, the complementary solution must also be antisynmetric . These 
complementary solutions are, respectively, 

v ) _ g cosh erg cosh grr) 

lk'cs cosh er cosh gr 



(v 


) = E 

lk ca 


sinh err] + ^ sinh grt} 
sinh er sinh gr 


(49) 


E - [sinh e r (l - f fc (£)d£ 


F 



[sinh g r(l - £]] 




(50) 


(51) 


The Vjjj. solutions can be written explicitly in real form (by proper 

combinations of particular and complementary solutions) and are given in 
equations (Bl) to (B6) in appendix B. Values of v^ for several R 

(or Ra) are presented in table II. 


Velocity and Temperature Distributions 

Now that the various v^ are known explicitly, they can be Inserted 
into equation (46) to yield Vj_. The sum of_ v 0 (as given by eq. (38)) 
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and v-^ then forms the solution which Includes frictional heating effects 

to a first order. By means of the various transformations made, the solu- 
tions of the given boundary-value problem (eqs. (22) to (25)) or the 
dimensionless velocity distributions to the zeroth and first approximation 
are, respectively, 


k a 

^0 " Ra 


-a 


v 00 + 


(m - l) Ra 1 / 2 C 


V 01 


(m + l) Ra 1 / 2 G 


v 02 


+ a 


(52) 


u 


1 Ra 


[64 (v 0 + v-j_) + ccK 

2n„n2 


+ V. 

- ^0 + 64Ra 
2 


(m - lV^RaC* . (m + l) 2 RaC : 

"t* a O *T* A 

4 


*12 + 4 V 13 -(“-!) “Ra' • -14 


a 2 v^Q-dRa d / 2 C(m + l) v -j -| 

.V^cv, 


(m 2 - l)RaC 2 
+ * 2 V 15 


(53) 


where the v Q j (j = 0, 1, 2) are given in equations (40 ) to (42) and the 
v-jjj. (k = 0, 1, 2, 3, 4, 5) are given in appendix B. 

In principle, higher-order approximations could be obtained by con- 
tinuing the procedure described. However, the results become very un- 
wieldy, Therefore, beyond the range of applicability of the zeroth- and 
first-order approximations (that is, in the range of large frictional 
heating effects), the complete boundary- value problem should be solved 
numerically^ some discussion of these numerical results relative to the 
zeroth- and first-order approximations will be presented subsequently. 


To determine the temperature distributions, recall that they are 
related to the velocity distributions by equation (17 ) and that 

„ Ra 1 / 2 , „ Ra 1 / 2 ^ 

V 00 T~ V 02 and v 02 = -~T~ 30 that 


'00 


T 0 - CK A - ug 



“172 av 02 + 2 ( m “ ^^Ol + 

Ra ' L 


(m + l)C 
2 



( 54 ) 
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and 


K A 

T = CK. - u" » T a - — — 
1 A 1 0 16Ra 


a ^ V l 0 " a * Ra ' 1 / 2 C(m + 1) v!^ 


+ - ^.^ggf V J 2 + (» * vls .(».!) oSal/2 Cv" 14 


1 (m - 1)W "I 

2 15 J 


(55) 


where the primes on the u-fnnctions denote differentiations with respect 
to y and on the v-funetions denote differentiations with respect to tj, 
and the explicit forms of the second derivatives appearing in equations 
(54) and (55) are presented in equations (B7) to (B13) of appendix B. 

Thus, solutions of the original boundary-value problem in terms of 
u and T are known to zeroth- and first-order approximations, that is, 
neglecting frictional heating and including its effects to a first order, 
respectively. These solutions are, however, given in terms of the param- 
eter C (recall that m is also a function of C through eq. (25)). 
Therefore, to relate C to the physical problem the dimensionless mass 
flow in the channel is defined as 


M 



Vy 


(56) 


Neglecting aerodynamic heating, equation (56) becomes (using the zeroth- 
order approximation as given by eq. (52)) 


M = 


, K A I" a/ 2 a (sin 2r + sinh 2r\ 

*0^ Ra j^l/ 4 \posh 2r + cos 2r / 


+ RaV 4 ^( m + l) ( sin 2r - sinh 2r 
2 V 


y/2 


cosh 2r + cos 2r 


) 


(57) 


or 


, 1 

(cosh 2r + cos 

2r\ 

MRa 

a j \fz a/si n 2r + sinh 2r\ 

= 2Ad f- n 1/4 

a/2 Ra 1 

\^in 2r - sinh 

2r/ 

_ k a 

l/4\cosh 2r + cos 2 rj 


( 58 ) 
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(For large Ea the expressions in parentheses in equation (58) reduce 
to 1 and -1, respectively.) Substitution of G as determined into equa- 
tions (52) to (55) then yields the velocity and temperature distributions 
for any case if M, ft^/Ad, Ra, cc,, and K A are known. Note, then, that 

the mass flow is an independent parameter of the problem, and hence could 
be due in part to a forced flow. For a given configuration and fluid the 
last four of these parameters are specified, so that equation (58) re- 
lates C to the mass flow M, that is, the temperature level and the 
mas s flow through the channel are related, as is only reasonable in flows 
of this type. 

Asymptotic solutions . - The solutions presented in the previous sec- 
tions are valid (to the proper order of approximation) for all values of 
the parameters of the boundary- value problem. However, from physical 
considerations, it can be seen that in many practical occurrences of the 
phenomenon under consideration, the parameter Ra may become very large 
(of the order of 10 4 and higher) . It is therefore appropriate to examine 
the asymptotic character of the boundary- value problem. To this end it 
is convenient to write equations (22) to (25) as 


JL_ 

Ra 


_iv 

u 


K A C 2 

rT* fa ') 2 + u 



(59) 


u(0) =* u(l) = 0 (60) 

u" (0) = 1 u" (l) = m (61) 

where 

u = CK^u (62) 

For very large Ra, equation (59) is of the boundary layer type (see ref. 
10) . Therefore the velocity and temperature profiles will have very large 
gradients near the walls, and thus the asymptotic solutions will yield the 
velocity and thermal boundary layers with essentially constant conditions 
given by the in viscid solution u^. = ct/CRa in the center of the channel 

associated with large Ra flows. Hence, expanding the coordinate normal 
to the wall, as Is done in boundary layer theory, requires that 

y = Ra 1 / 4 y (63) 


u = Ra" 1//2 V 


(64) 
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The asymptotic forms of equations (59) to (6l) are 


CK 


yyyy 


+ v -sr c V -~V = 0 

Ka CRa ' 


v(0) 


_ 1 for the wall at temperature T w _ 

V — (0) = m = { ° 

TV 1 m for the other wall at temperature 


V(») = V^(CO) = 0 


(65) 

( 66 ) 

(67) 

( 68 ) 


where the subscripts denote differentiation. Note that the conditions 
expressed by equation (68) replace the boundary conditions at the second 
wall and require that influences of one wall do not affect the other. 

Thus y can be considered as the coordinate normal to the first (or left- 
hand) wall (that is, the one corresponding to Y = 0) and -y will be 
the coordinate normal to the other wall, which is at Y =* d. Hence, use 
should be made of the proper part of equation (67) in each solution. 

From equation (65) it can be seen that the frictional heating effects 
will be negligible for large Ra unless CK^/Ra (or Cy3f x d/cp) iB at 

least of unit order of magnitude. It should perhaps be pointed out here, 
in contradistinction to the qualitative discussions In references 1 and 
12, that the frictional heating is important only If essentially the 
ratio of K (based on any appropriate temperature) to Ra is of unit 
order of magnitude or larger, as can, in fact, be verified in general. 

Thus it should be noted that the discussions in those references hold 
specifically only if Ra is of unit order or smaller or If Ra does not 
appear as an explicit parameter of the problem (as in ref. 1, for example). 
For the range of conditions and physical properties of fluids being con- 
sidered, It is unlikely that CK^/Ra will be of unit order for a flow 

generated In a gravitational field alone with large Ra. Therefore, un- 
less the natural-convection flow is being generated by a body force con- 
siderably stronger than gravity, the ratio of the volumetric expansion 
coefficient to the specific heat at constant pressure is unusually large, 
or there is considerable forced flow (to increase C), the frictional 
heating effects will not be important for large Ra. (Of course, the 
possibility always exists that some unusual fluid will be employed whose 
physical properties are such that CK^/Ra will be of such a magnitude 

that the frictional heating effects will be important for large Ra even 
in a gravitational field. Liquids near their critical state may be rep- 
resentative in this respect (see ref. 13)). Furthermore, it can be seen 
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from the same equation that the effects of the heat sources will "be im- 
portant only if a/CRaV^ is of unit order of magnitude for large Ra. 
This is physically reasonable. 


Since the frictional heating effects are negligible for large Ra, 
unless CK^ is very large these effects will be neglected in this sec- 
tion. (A method of successive approximations similar to that described 
in the previous sections could be applied without difficulty to eg. (65) 
if these effects are of consequence.) Therefore, letting 

V = V - a/CRa 1 / 2 , equations (65) to (68) become, respectively. 


V + V 

yyyy 


= o 


(69) 


V(0) = - 


a 


CRa' 

*yy(°> = “ 


172 


(70) 

(71) 


VH = V- = 0 


(72) 


The solution satisfying equations (69) to (72) is 


sin 


To find the temperature distribution the second derivative of equation 
(73) is necessary, and this is given by 



a 


-r cos 

V2 CRa 1 / 2 V*/ 


_2_ 


(73) 



-y/Vs 

e 





( 74 ) 


For large values of the Rayleigh number Ra, then, the dimensionless 
velocity and temperature distributions are given by 


and 


u 


a 





(75) 


T a = GK A C 1 - 


(76) 
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Solution of Special Case Simulating an Enclosed Channel 

There is considerable interest in the natural-convection flow in a 
completely enclosed rectangular region. In reference 7 , the natural- 
convection flow of fluids containing heat sources between two parallel 
planes is considered (els is the case here also); but there, in order to 
simulate an enclosed rectangular region, it is further specified that the 
net mans flow be zero and that the walls be at the same temperature. This 
problem is treated in an approximate manner in reference 7 in that the 
velocity distribution was postulated without regard to the equations of 
motion and, hence, it would be desirable for comparison purposes to obtain 
a more exact solution for this special case from the solutions found 
herein. 

For zero net mass flow in the channel the parameter C can be deter- 
mined directly from the other parameters of the problem and equation (58) . 
Thus 



a 

a/2 Ra 1 / 4 


/cosh 2r + cos 2r\ 

— - 

rj% /sin 2r + sinh 2r\ ^ 

\sin 2r - sinh 2 rj 

Ra l/4 \cosh 2r + cos 2r / 


(77a) 


where the superscript denotes the zero net mass flow and the subscript 
designates the order of the approximation. It is interesting to note 
from this equation that for no internal heat sources (a =« 0), zero net 
mass flow in the channel can be obtained when C ** 0 W ^/ 2Ad. However, to 

obtain the solution for the special case simulating flow with heat sources 
in a completely enclosed region with walls of equal temperatures from the 
solutions presented in the previous sections, must be zero in equa- 

tion (77a), and therefore 


C 


0 

0 


= a 1 

(cosh 2r + cos 2r\ 

\/2 . j 

/sin 2r + sinh 2r ^ 


a/2 Ra 1 / 4 

(sin 2r - sinh 2r / 

La 1 /* ' 

(cosh 2r + cos 2r J 


(77b) 


The velocity and temperature distributions for this special case are then 


obtained by replacing C in equations (52) and (54) by 


as given by 


equation (77b). Further physical significance of the zero net mass flow 
case can be inferred from equation (13) . Integration of this equation 
over the channel cross section shows that for no net mass flow all the 
heat generated internally (by heat sources and by aerodynamic heating if 
the latter is significant) in a given cross section must be transferred 
to the walls. 
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Computations made for this case including frictional heating yielded 
no appreciable deviations from the zeroth-order results. 


Solutions for Case of Convective Inversion 

An interesting aspect associated with the natural-convection 
phenomenon is that the Grashof number can change sign; this implies a 
reversal in the flow direction and is referred to as convective inver- 
sion (see p. 109, ref. 8). The sign of the modified Grashof number in 
the present paper can be changed in one of three ways: (l) by a change 

of the sign of the longitudinal temperature gradient cffi/dX = A, (2) by 
a change in the direction of the generating body force, and (3) by a 
change in the sign of the volumetric expansion coefficient as occurs 
near the critical state of a liquid (see ref. 13). 

Since the modified Grashof number appears in the parameters 

and Ra connected with the problem considered herein, the effects of 
sign changes of the Grashof number in the solutions should be studied. 
From its definition, is proportional to A, f x , and (3^ so that 

only the first will alter its sign. Note further, however, that the 
modified Rayleigh number Ra is essentially the product of the Prandtl 
number and a modified Grashof number, that is, a Grashof number which 
depends on the product (Ad) of the longitudinal temperature gradient and 
the distance between the plates. Hence, any one of (l), (2), or (3) 
given in the preceding paragraph will lead to a change in the sign of 
Ra. For negative Ra the solutions as given in the previous sections 
do not apply, and hence the foregoing boundary-value problem (eqs. (22) 
to (25)) would have to be solved with negative Ra. These solutions 
can be readily obtained, but it is found that with frictional heating 
neglected these solutions change character with changes in Ra and that 
there exist critical negative values of Ra for which the solutions 
become meaningless. In an attempt to explain these unusual results, 
further interpretation of the problem must be made. Reexamination of 
the meaning of negative Ra shows that not only changes in the body 
force direction and sign of the volumetric expansion coefficient but 
also a change in the sign of the longitudinal temperature gradient A 
can lead to negative Ra. If the negative Ra is attributed to the 
last cause, the physical interpretation of the unusual mathematical 
results pointed out becomes clearer, because a negative A implies that 
the fluid is being heated from below and this situation leads to a 
"Rayleigh-type" instability of natural-convection flows due to the 
"p ilin g of heavy fluid on lighter fluid." Analogous interpretations, of 
course, also follow directly for changes in the body force direction and 
in the sign of the volumetric expansion coefficients. Natural -convection 
flows heated from below between horizontal plates have been studied ex- 
perimentally in some detail (see refs. 8, 14, and 15, for example), and 
it was found that the flow does indeed change character (into cellular 
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motion) for certain critical values of the Rayleigh number. Hence, it 
is believed that the critical values of Ra found from the linearized 
analysis (in that the aerodynamic heating was neglected) for negative 
Ra may be analogous to those observed in actual cases. However, since 
this instability leads to these additional complications and should be 
further investigated (perhaps using the true nonlinear eqs.), the case 
of negative Ra will not be treated further herein. It should, how- 
ever, be kept in mind that if in an actual case of the configuration 
considered herein the Ra is negative, the flow and heat transfer will 
not be as predicted in this paper but should be expected to exhibit a 
behavior pertinent to the "unstable-type" flows. 


RESULTS AN D DISCUSSION 

Velocity and Temperature Distributions 

The relations between the actual and dimensionless velocities and 
temperatures as determined from the various transformations in the 
analysis (see eqs. (16)) are 



where U and Q denote the actual and u and T, the dimensionless 
quantities. For a given fluid, configuration, heat-source Intensity, 
and mass flow, the velocity and temperature distributions can be com- 
puted from equations (52) to (55) (for Ra > 0); and for zero net mass 
flow and the walls of the same temperature, by applying equation (77b) 
to equations (52) and (55). These computations will be accurate within 
the limits of the method of solution; that is, for moderate and small 
values of Ra the solutions yield results of reasonable accuracy for 
small CK^, and for large Ra the zeroth-order approximations or, even 

more simply, the asymptotic solutions will give answers valid for all 
CK^. The range of applicability of the various solutions presented 

herein will be discussed more fully subsequently. 


Because the solutions were obtained in the convenient forms (eqs. 
(52) and (54), for example) wherefrom the qualitative effects of the 
various parameters associated with the problem can be studied, and since 
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tabular values of the universal functions are presented to facilitate 
computations for any specific case , no extensive detailed calculations 
will he given covering the entire range of values taken on hy the param- 
eters. Representative velocity and temperature profiles were, however, 
calculated for Kj^ * 10, m = -1, 1, and 2, a, =* 0, 10, and 100, and 
Ra = 10, 10 2 , 1600, and 10 4 . In addition, the parameter C was given 
the value -1 in a3 3 the computations except those for the case simu- 
lating flow in a completely enclosed region. In this way the relative 
influences of the other parameters are just as apparent, hut the number 
of computations is greatly reduced. The results of these computations 
are presented in figures 2 to 9. The contents of each specific figure 
(numbers 2 to 9) are listed in the following table: 


Ka 

= 10 


^0* 

a l 



T 0> 

T 1 


m 

a 

Ra = 10 

10 2 

1600 

10 4 

10 

10 2 

1600 

10 4 

-1 

0 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 


10 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 


100 

2(a) 

3(a) 

4(a) 

5(a) 

6(a) 

7(a) 

8(a) 

9(a) 

1 

0 

2(h) 

3(h) 

a 4(h) 

5(h) 

6(h) 

7(h) 

a 8(h) 

9(h) 


10 

a ^2(h) 

3(h) 

4(h) 

5(h) 

a ^6(h) 

7(h) 

8(h) 

9(h) 


100 

2(h) 

3(h) 

4(h) 

5(h) 

6(h) 

7(h) 

8(h) 

9(h) 

2 

0 

a ^2(c) 

3(c) 

4(c) 

°5(c) 

a ^6(c) 

7(c) 

8(c) 

C 9(c) 


10 

2(c) 

3(c) 

4(c) 

5(c) 

6(c) 

7(c) 

8(c) 

9(c) 


100 

2(c) 

3(c) 

4(c) 

°5(=) 

6(c) 

7(c) 

8(c) 

C 9(c) 


a Includes results for u or T. 
^Includes results for u( 2 ) or t( 2 ). 
c Includes results for u^ or T a . 


For each triplet of parametric values (m, a, Ra) the profiles were com- 
puted with frictional heating neglected (hy eqs. (52) and (54) and de- 
noted hy Uq and T Q on the figs.), with frictional heating included 

to a first approximation (hy eqs. (53) and (55) and denoted hy u-j_ 

and Tq), and in several specific cases with frictional heating completely 

accounted for (hy numerical solution of eqs. (22) to (25) using a Card- 
Programmed Electronic Calculator and denoted hy u and T ) . For 




22 


KACA TN 3141 


Ra = 10 4 the asymptotic solutions (given "by eqs. (75) and (76) and 
denoted "by u a and T a ) are also included for =* 10 , m =* 2, and 
a = 0 and 100. The asymptotic solutions were computed from each in- 
dividual wall to the channel center and then faired in to join smoothly. 
From these computations any qualitative trends obtained by examination 
of the solutions can be further substantiated and, in addition, some 
definition of the range in which the frictional or aerodynamic heating 
exerts a large influence can be made. Calculations were also made for 
the special case simulating a completely enclosed region in which there 
is no net mass flow and the walls are at the same temperature by applying 
equation (77) to the appropriate solutions. These curves are given in 
figures 10 and 11. The velocity and temperature profiles (particularly 
for Ra < 104) are qualitatively similar to those determined experi- 
mentally in reference 16. 

Effect of different wall -temperature configurations (m varying) 
and heat sources (a varying) . - From equations (52) to (55) and their 
related ■universal functions, it can be seen, as expected, that an in- 
crease in the wall temperature parameter m or an increase in the heat- 
source parameter a results in larger velocities and higher tempera- 
tures. These trends together with that of increasing net mass flow, as 
represented by the area under the u-curves, with m and a can be ob- 
served on figures 2 to 9. It can also be seen from figures 6 to 9 that 
if sufficient heat is generated by the heat sources, the direction of 
heat transfer will be changed. In agreement with the statements made in 
the section dealing with the zero net mass-flow case, note from fig- 
ures 2(a), 3 fa), 4(a), and 5(a) that if aerodynamic heating is neglected 
for m •=* -1 (since C = -l) and a » 0, there Is no net mass flow. In 
general, the velocity distributions become more symmetrical with the 
larger a (see figs. 2(a), 3(a), 4(a), and 5(a), for example) because 
the heat added uniformly by the heat sources counteracts any asymmetry 
imposed by the wall thermal conditions. 

Effect of the modified Rayleigh number (Ra) . - Examination of 
the solutions (eqs. (52) to (55)) shows that the' velocities and tem- 
peratures decrease with increasing values of the modified Rayleigh 
number Ra. This trend can also be seen by comparison of corresponding 
curves in figures 2 to 9, and even by comparing with the curves in refer- 
ence 1 which are for Ra « 0. For large Ra it can be seen from fig- 
ures 5(b), 5(c), 9(b), and 9(c) that the velocity and temperature pro- 
files take on a "boundary-layer form." Asymptotic solutions computed 
for Ra = 10 4 , = 10, m = 2, and a = 0 and 100 are also presented 

on figures 5(c) and 9(c), and these very closely approximate the more 
exact solutions. Hence, for large Ra the asymptotic solutions can be 
employed to yield reasonable results much more simply. For the case 
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where m = -1 and a = 0, increasing the modified Rayleigh number changes 
the temperature distribution from essentially the conduction profile 
(that is, an almost linear distribution) at Ra = 10; the increased 
effect of the convection is then apparent for the larger values of the 
Rayleigh number. It is interesting to note that for = 10, a =■ 10, 
and m = -1, 1, and 2, changes in the Rayleigh number can so affect the 
temperature distributions that the heat flow direction from one or both 
the walls can be altered. (Compare corresponding parts of figs. 6 
to 9.) This point will be more graphically portrayed in the subsequent 
discussion of Nusselt numbers. 

The velocity and temperature distributions (see figs. 10 and ll) 
for the special case considered herein of zero net mass flow and walls 
at the same temperature are not in general appreciably altered in shape 
by increases in Ra although the velocity peaks vary inversely with the 
Rayleigh number. The shape of the velocity and temperature profiles is 
seen to be qualitatively the same as that assumed in reference 7. 

Effect of frictional heating . - By comparing the profiles presented 
in figures 2 to 9 computed by neglecting frictional heating (denoted by 
the subscript zero) with those computed including the aerodynamic heat- 
ing to a first order (denoted by the subscript unity), the effect of the 
aerodynamic heating on the velocities and temperatures can be studied, 
numerical solutions obtained of the complete boundary-value problem 
(eqs. (22) to (25)) in which the frictional heating was entirely taken 
into account are also included (with no subscripts) (see preceding table) 
on figures 2(b), 2(c), 4(b), 6(b), 6(c), 8(b), 10, and 11 for comparisons 
with the approximate solutions. 

In accord with the discussion on the asymptotic solutions, it can 
be seen that when K^_ is small compared with Ra (recall that C = -1 
in these calculations) the aerodynamic heating effects are negligible. 
Since no computations were made herein for > 10, the computations 
made for Ra = 1600 and 10^= show no deviation between the zeroth and 
first approximations, and these are also coincident with the numerical 
solution (see figs. 4(b) and 8(b)). Hence, in the range K^/Ra « 1, 
the zeroth-order approximations will yield accurate results; if, in 
addition, Ra>> 1, the asymptotic solutions provide a simple means of 
obtaining the velocity and temperature profiles. Note that the param- 
eter Ka serves merely as a scale factor in the zeroth-order solutions. 

In the range where and the Rayleigh number are of the same 

order of magnitude, the frictional heating affects the results to greater 
or lesser degree depending on the particular amount of heat addition as 
specified essentially by the parameters m and a (figs. 2, 3, 6, 
and 7) and, hence, the first-order approximations should be employed in 
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this range. Comparison, of the first-order solutions with the several 
numerical solutions of the entire "boundary-value problem shows close 
agreement for the cases computed. Hence, unless the conditions are 
more severe than the most extreme conditions for the range of parameters 
considered herein as represented "by K^/Ra = 1, m = 2, and a = 100, 
the first-order solutions will yield results which include the effect 
of aerodynamic heating of reasonable accuracy. In this range where its 
effects are important, this frictional heat, of course, acts just as do 
the heat sources and leads to increased velocities and altered tempera- 
ture profiles and, consequently, different heat-transfer rates (figs. 2, 
3, 6, $nd 7)v. 

For the special case simulating flow in a completely enclosed re- 
gion, the aerodynamic heating did not affect the results appreciably over 
the range of parameters under consideration. A numerical solution com- 
pletely including the effects of aerodynamic heating was made for 
Ra =* 10; it can be seen in figures 10 and 11 that this solution coin- 
cides with the zeroth-order solution, which neglects the effect of fric- 
tional heating. 

The complete consideration of frictional heating (as in the numeri- 
cal solutions) for the problem discussed herein, just as for the case 
(essentially Ra » 0) reported in reference 1, leads to the two results 
(a) that there exists a critical set of conditions beyond which no solu- 
tions exist, and (b) that where solutions exist there are two solutions 
for every set of admissible parametric values. Examples of these second 
solutions are presented in figures 2(b), 2(c), 6(b), and 6(c) (denoted 
by u(2) and t( 2)), and it can be seen that the velocities are more 
than 10 times as large as the first solutions and the temperatures are 
much greater than the corresponding first solution temperatures. These 
last unusual results cannot be predicted from the solutions obtained by 
successive approximations as described herein, but are found from numeri- 
cal solutions obtained by means of a Card-Programmed Electronic Calcu- 
lator. At present the significance of the second solutions is not ex- 
plained, although it is felt that they are intimately connected with the 
unique regenerative action of the frictional heating in natural convec- 
tion. The existence of the critical conditions appears to be similar 
to the thermal choking phenomenon. 


BEAT TRANSFER 
Nusselt Numbers 

The heat-transfer coefficients for the natural-convection phenomenon 
treated here can be expressed in terms of Nusselt numbers. The Nusselt 
number is here defined aB 
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where the double subscript signifies that the temperature gradient is 
to be evaluated at either Y = 0 or Y = d, depending on the wall under 
consideration. 

In terms of the dimensionless quantities, 



(80) 


The temperature gradient can be found from the zeroth-order solutions, 
and the Nusselt n umb ers can be computed on this basis from 


Nuoo = 


R 


jlA 


2 2 

cosh r + cos r - 1 


a(sinh 2r + sin 2r) 


+ c ( IlH - 1 ) Ra1 ^ 2 (siuh 2r - sin 2r) + — C - ( m r 1 . ) Ra L — _ — (sinh 2r + sin 2rl 

J 2(cos r - cosh r) J 


(81) 


_ 1 J -1 

NU 01 = e- 1/4 ] ■ 2 2 “ 

2v2 Ra ' t_cosh r + cos r - 1 |_ 


L(sinh 2r + sin 2r) 


+ (sinh 2r - sin 2r) 


C(m-l)Ra 1 / 2 
2(cos 2 r - cosh 2 r) 


(sinh 2r + sin 2ry 


J 


(82) 


where the first subscript denotes that zeroth-order approximation is 
used, and the second denotes the wall with which the Nusselt number is 
associated. (Eqs. (81) and (82) are, of course, specifically for Ra 
positive.) Note that these zeroth-order Nusselt numbers are independent 
of K^. When CK^ is of the same order of magnitude as Ra, the zeroth- 

order approximation has been shown to be inaccurate; therefore improved 
Nusselt numbers can be obtained by using the appropriate t (first- 
order approximation) solutions and, for the same conditions as for 
equations (8l) and (82), the Nusselt numbers can be computed from 


Nuio = Nuqo + ~ wtt P&L " CoRal/2 ( m + ijQg + ^ 7 

K4»fo Ra'/ L_ 


64^" Ra 

n 2 /■„ , ,\2 


l)^n t + 0(m . 1)aai V^ - gflsll 11 EafJg" 


+ 


4 


2 


(83) 
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Nun i =. Hu™ — — =-7— - CaRa^/^fm 

11 ^ 641/2 Ra 7 / 4 L 

+ c2 (*^ )2Ba a t + C(m - 

where the various 2^ in equations (83) and 
(B43) to (B48) in appendix B. 


♦ ift 2 4 

+ c?& 2 - . M (84) 

(84) are given by equations 


Computations of the Nusselt numbers were made over a range of 
values of Ra from equations (8l) to (84) for = 10, m = -1, 1, and 
2, and a =* 0 and 10 where again C was taken to he -1, and the re- 
sults Eire presented in figures 12 and 13 for the wall at Y =* 0 and 
Y = d, respectively. The figures show that for the wall at Y => 0 the 
Nusselt numbers decrease, in general, with increasing Ra, and for the 
wall at Y =» d the Nusselt numbers increase with the modified Rayleigh 
number except for the case where m = -1. For = 10 and a 10, 
the Nusselt number changes sign with increases in Ra because of the 
variation of the temperature profile with Ra, as was previously noted. 

Figures 12 and 13 also demonstrate clearly the effect of the aero- 
dynamic or frictional heating. This effect, in accordance with all that 
preceded, is extremely pronounced for low values of the Rayleigh number 
(that is, when K/Ra is of unit order of magnitude). 

Flow in an enclosed region . - For the special case simulating flow 
in a completely enclosed region (M = 0, m = 1, and a ^ 0), the calcula- 
tions for the temperature profile were extended over a larger range of 
Ra and plotted in figure 14 as the ratio $ of the temperature differ- 
ence to that for pure conduction, as was done in reference 7, The con- 
duction temperature difference used in $> is the channel center-to- 
wall difference subject to uniform heat generation by sources and is 
equal to aKpJ 8. The temperature profiles computed in an approximate 
manner in reference 7 Eire compared with those computed more exactly by 
the method reported in this paper, and it can be seen that the dis- 
crepancy becomes quite apparent for values of Ra = 10^ and above. 

The Nj in reference 7 is related to Ra by Ra = ^ Njj hence, for 

most liquids Nj and Ra are identical. If T / 1, T is merely a 
scale factor. It can be seen from figure 14 that for Ra < 10^, the 
temperature gradients at the wall are all almost identical, and hence, 
even if the temperature profiles themselves were not identical, the 
heat transfer computed by the two methods would be in reasonable agree- 
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ment. The variations of the dimensionless temperature variable 
(ratio of eenter-to-wall temperature difference to that for pure con- 
duction) as used in reference 7 with Ra as given by the two methods 
are also compared in figure 15, and hence the quantitative limits of 
the approximate method can be seen. 


CONCLUDING REMARKS 

An analysis was made of the flow subject to body forces between two 
parallel plane surfaces oriented in the direction of the generating body 
force along which the temperature is specified to vary linearly. The 
solutions for natural convection and those for combined natural and 
forced convection were found to be essentially the same. It was found 
that a modified Rayleigh number (product of Prandtl and modified Grashof 
numbers) in addition to the parameter K^ was of significance in this 
problem. Solutions for the velocity and temperature distributions are 
given in terms of "universal" tabulated functions. Detailed velocity 
and temperature profiles were computed and it was found that, in general, 
the velocity and temperature differences increase with the wall tempera- 
ture parameter and with additional heat due to heat sources . The veloc- 
ities and temperatures decrease with increasing values of the modified 
Rayleigh number. When the ratio of CKfl. to the modified Rayleigh num- 
ber is of unit order of magnitude, the frictional or aerodynamic heating 
appreciably affecus the velocity and temperature distributions. Asymptotic 
solutions for large Ra are presented which make computations in this 
range relatively simple. For any given set of the parameters, complete 
consideration of frictional heating implied the existence of two flow and 
heat-transfer states and implied that no solution exists beyond certain 
critical values. 

Consideration was given to a special case simulating the natural- 
convection flow of fluids with heat sources in a completely enclosed 
region with the walls at the same temperature. Computations from the 
solutions for this special case demonstrated that an approximate method 
developed in another paper should yield reasonably accurate results as 
long as the modified Rayleigh number is less than 10^. The effects of 
aerodynamic heating were found to be negligible for this case. 

Study of the convective Inversion aspect of the present problem led 
to the Inference that for negative values of the Rayleigh number, addi- 
tional complications arise because of an instability (due to heating from 
below) of the flow which must be more thoroughly investigated. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 29, 1953 
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APPENDIX A 
SYMBOLS 

The following notation is used in this report: 

A longitudinal temperature gradient 

a.^ constants in eq. (39); i = 1, 2, 3, 4 

B ± constants defined by eqs. (B14) to (B35)j i = 1, 2, 3,. ..22 

C constant in eq. (17) 

constants defined by eqs. (B36) ; i = 1, 2, 3, ...6 

C° constant defined by eq. (77b) 

C constant In eq. (14b) 

Cp specific heat at constant pressure 

c v specific heat at constant volume 

D^ constants defined by eqs. (B37); i = 1, 2, 3, ...10 

d characteristic length (specifically distance between plates) 

E constant defined by eq. (50) 

E^ constants defined by eqs. (B38) ; i ^ -1, 2, 3, ...6 

e constant, (i+l) 

F constant defined by eq. (5l) 

f-jj- negative of X-component of body force per unit mass 
» 

G Green's function 

4/2 

Gr^ modified Grashof number , pf^ Ad /y 
g constant, (i-l) 
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h. 

Hi 

J 

K A 

k 

Hi 

M 

Mi 

m 

m 

Ru 

% 

P 

Pr 

Q 

Qi 

R 

Ra 

r 

T, T* 
U 
u 
u 


heat -transfer coefficient 

constants defined by eg.. (B39); i = 1, 2, 3, 4 
constant, Ra 1 /^CK^/64 

dimensionless parameter . Pr Gr* ~ — 

Cp 

thermal-conductivity coefficient 
constants defined by eqs. (B40)j i = 1, 2 


dimensionless mass flow 

constants defined by eqs. (B4l ) ; i = 1, 2, 3, 4 
constant defined by eq. (25) 
constant defined by eq. (64) 

Nusselt number, hd/k 

modified Rayleigh criterion as given in ref. 6, yRa 
pressure 

Prandtl number, CpP./k 


heat due to heat sources 

constants defined by eqs. (B42)j i = 1, 2, 3,... 8 
constant, Ra/l6 

modified Rayleigh number, — Pr Gr^_ 


constant, r i/4 /V5 

temperature 

velocity 

dimensionless velocity 
dimensionless velocity, u/CK A 



30 


NACA TN 3141 


V 

V 


v 

V 

X 

Y 

y 

y 

a 


P 


r 

n 

e 


A 

M- 

5 


P 

T 

$ 



dimensionless velocity defined by eg. (6l) 


dimensionless velocity, V - 


a 


CRa 1 / 2 

dimensionless velocity, v/64 
dimensionless velocity defined "by eg. (26) 
longitudinal coordinate 
transverse coordinate 
dimensionless transverse coordinate 
dimensionless transverse coordinate, EaV*y 
dimensionless heat-source parameter, Qd/kA 


coefficient of volumetric expansion, p 


[WpI 


<5r 


J p 


ratio of specific heats 
dimensionless coordinate, 2y - 1 
temperature difference, T - T VQ 

constant , -a KpJ 64 

absolute viscosity coefficient 
kinematic viscosity coefficient 
dummy variable 
density 

dimensionless temperature difference 
d ime nsionless temperature difference, 

dimensionless center-to-vall temperature difference, 
constants in egs. (80) and (8l); i = 1, 2, 3, ...6 
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Subscripts: 

a asymptotic solution 

ca complementary antisymmetric solution 

cb complementary symmetric solution 

I inviscid solution 

n order of approximation 

P particular solution 

vq conditions at y = 0 

w-j^ conditions at y = 1 

0 zeroth-order approximation 

1 first-order approximation 
Superscript : 

(2) second flow and heat -transfer state 
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APPENDIX B 


FIRST -CKDER SOLUTIONS 

The explicit forms of the first-order solutions to "be used in equa- 
tion (46) are: 


*10 = — ^70^1 cosh rq cos rq + Bg sinh rq sin rq + B3 cosh 2 ttj cos 2rTj 

* L 

+ B4 sinh 2rr) sin 2^1 + 65(008 2tt] - cosh 2rq) + 15 B^J (Bl) 

Vn = — i—fBg cosh rt] cos ri] + By sinh rq sin rq - B4. cosh 2rr\ cos 2rq 
Ra 3 / 2L 

+ B3 sinh 2 tt) sin 2 tt] - 15 B^J (B2) 

*12 - — iy r 2’[j B 8 cosil rT l cos rT l + b 9 rT l si* 1 rT l 

+ B 10 cosh 2 rq cos 2 rr] + B-^ sinh 2 rq sin 2rTi 

+ B 12 (cosh 2rTj + cos 2rTi) - 15 B-^JJ (B3) 

^jg|Bi3 cosh rq cos rq + B^ sinh ri} sin rrj 


Ra" 


- B 3 cosh 2rq cos 2rr) - B^ sinh 2rT} sin 2rT] 
+ B 5 (cos 2rr\ - cosh 2rq) - 15 B^j 


v 14 - 




B^. sinh rq cos rT] + B n K cosh rT] sin rq 


16 


+ B 17 sinh 2rr] cos 2rr| + B 1Q cosh 2ri] sin 2rT] 
+ B ig sin 2rt] + B 2Q sinh 2rrj| 


(B4) 


(B5) 
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v, r = — \ t ? [b o , sinh rq cos rq + B co cosh rq sin rq 
15 Ra°/ ^ l_ 21 “ 

- B^_g sinh 2rT) cos 2rT] + B-jj cosh 2rq sin 2rq 
+ B 2 q sin 2rq - B^g sinh 2rqJ (B6) 

The constants B^(i = 1, 2, ... 22) appearing in the preceding equations 

are readily computed for a given Ra. In the subsequent section of this 
appendix these constants are written explicitly in a form suitable for 
reasonably rapid computation. 

The second derivatives appearing in equations (54) and (55) are 


-Ra 1 / 2 


[cosh r sin r cosh rq sin rq 


01 4(cos 2 r - cosher) u 

+ sinh r cos r sinh rq cos rqj (B7) 

v 10 = c ° 8il rT l cos rT l - Bp sinh rq sin rq 

+ 4B^ cosh 2rq cos 2rq - 4Bg sinh 2rq sin 2rq 
- 2Bg(cos 2rq + cosh 2rq/j (B8) 

v^i = cosh rq cos rq - Bg sinh rq sin rq 

+ 4Bg cosh 2rq cos 2rq + 4B^ sinh 2rq sin 2rqJ (B9) 

v 12 = ^a[^9 cos ^- rI l cos rT l " Bg r B sin rT l 


+ 4B-[ -] cosh 2rq cos 2rq - 4 B^q sinh 2rq sin 2rq 
+ 2B^g(cosh 2rq - cos 2rq^ 

V 13 = 4^ [^14 cosl rT l cos “ ®i3 sinh rq sin rq 

- 4B^_ cosh 2rq cos 2rq + 4Bg sinh 2rq sin 2rq 

- 2Bg(cos 2rq + cosh 2rq)j 


(BIO) 
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v^ 4 = 4^[- b 15 cosh r1 l sin rT 5 + b 16 sirf:i r7 l cos rT| 

- 4B 1? cosh 2rrj sin 2rT] + 4 B-j_q sinh 2rrj cos 2rTj 

- 2B^g sin 2rrj + 2B 20 sinh 2rr]|j (B12) 

v 15 = 4^[_“ B 21 COSlL rT1 Sin rT1 + 522 SlD]:1 rT * C0S rT| 

+ 4Bi8 cosh 2rrj sin 2rr) + 4Br7 sinh 2rt] cos 2 tt) 

- 2B 2 q sin 2 tt \ - 2B ig sinh 2rrjj (B13) 

Constants for First-Order Approximation 

The constants appearing in equations (Bl) to (B13) are written in 
an expeditious form for confuting as follows: 


B n = - — 
1 % 


[(-3D 1 +D 2 )E 1 + (D 3 -3D 4 )E 2 + 2(D 3 +D 4 )E 3 + 2(0^^ 

+ ( b 3 +D 4)% “ ( b 1+ b 2)E 6 ]+ W 

Bg = - A[(D 3 -3D 4 )E 1 + (3D 1 -D 2 )E 2 + 2(D 1 +D 2 )E 3 - 2(D 3 +D 4 )E 4 


h. 

- (D 1+ D 2 )E 5 - (D 3+ D 4 )E 6 ]- (2hl ■ + -— 


Mi 


_ 32 h l h 2 

B 3 “ " 15 ~M^" 

16 &i* " ^ 
B 4 - JE ' 


Mi 


B 5= - 5 


16 ( h l 2 + t 2 2 ) 


Mi 


B 6 = ^E 3D 1- D 2 )E 2 + + <W E 6 + < D 1 +D 2 )E S 

(i-2 2 - - Hi 2 )] 


256 

15Mi 


(B14) 


(B15) 

(B16) 

(B17) 

(BIS) 


(B19) 
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128 h-, he 


B? = ^-^(-D 3 +3D4)E 2 + (5D 1 -D 2 )E 1 - (D-|_+D 2 )Eg + (%+D 4 )E^| + 

B 8 = - + (ZDi«7-®8 )C 2 - 2 ( D 5 +D 6) C 3 " 2 ( D 7 +D 8 )C 4 

448 h h 2 32 (M 4M ) 

+ (■ d 5“ I) 6 +2I ^4^5 + (-2Di+D7-D8)C 6 j+ 

Bg = - ^|^(2D-j_+Dy-Dg)C2_ + £-D 3 +Dg+2D 4 )C 2 - 2(Dy+Dg)C 3 + 2(Dg4dDg)C 4 

64(M 5 -M4.) 

+ ( -2D 1 +D 7 -D 8 )C 5 + (-D 5 +Dg-2D 4 )CgJ+ — 15Ms 


(B20) 


(B21) 


®10 = 


B H = 


B i2 “ at 


32 h-jh 2 
lSMg 

16(M 5 -M 4 ) 

15Mg 

16(M 5 +M 4 ) 


(B22) 

(B23) 

(B24) 

(B25) 


B^ = jj - “3D-^+D 2 )Ei + (D 3 -3D 4 )E 2 - 2(D 3 +D 4 )E 3 - 2(Dj4D 2 )E 4 

512 h-|_h 2 


+ ( d 3*®4) e 5 “ ( D l +D 2^ E 63 “ 


15M-, 


(B26) 


®14 = ^“[( D 3“ 3 ®4)El + (3D2~D 2 )E 2 - 2 (D-j_-®2 ) E 3 + 2(D3+D 4 )E 4 - (D 1 +D 2 )E s 


(D 3+ D 4 )E^] - [ hl 2 + Eh/) 


(B27) 


B is = ^ [( 2 Di+D ? -D q ) Qq. + (D 5 -Dg-2D 4 )Q 2 - (Dg-ffig)^ + (Dy+Dg)^ 

- (Dg sinh 2r)Q 5 - (D 10 sinh 2r)Qg - (D g sin 2r)Qy - (D 10 sin 2r ) t 3ej 
+ ^p[-20h 3 + 12h 4 + 18 sinh 2r - 6 sin 2rj (B28) 
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s 16 = ^■[(“ :d 5‘®6 +2D 4^1 + (2Di+ d 7“I ) 8)Q2 + ( D 7"®8^% + ( D 5"®6^^4 
2 

- (D^q sinh 2r)Qg + (Dg sink 2r)Qg - (D-^q sin 2r)Q 7 
+ (Dg sin 2r)Qgj+ ^=^ |l2hg - 20h 4 - 6 sinh 2r + 18 sin 2rj 


16h„ 


B 


'17 " 15M, 


B 


16h, 


18 = 


B 19 " ' 


b 20 = 


15M 3 

8 sinh 2r 
5M 3 

8 sin 2r 

®S 


B 21 " ^[(-SD^+Dq)^ + (D 5 -D 6 -2D 4 )Q 1 +_(D 5 +D 6 )Q 4 + (D 7 ^ 8 )Q 3 

u 

+ (Dg sinh 2r)Q7 + (D^q sinh. 2 t)Qq - (Dg sin 2r)Q5 
- (D^q sin 2r)Q^|+ rj-^-[2Qh 4 + 12hg + 18 sin 2r + 6 sinh 2rj 


B 22 = 4 t (l>5 ' I ’ 6 ' 2D4)Q2 + ( 2 D 1 + VV «1 - + < D 5-®6>«3 

+ (D 10 sinh 2r)Qy - (Dg sinh 2r)Qg - (D 10 sin 2r)Q 5 
+ (Dg sin 2r)QgJ- jj=^ jl2h 4 + 20h3 + 6 sin 2r + 18 sinh 2rJ 


(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 


(B35) 
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where 


2 4 

C-j_ = - j sinh 2r sin 2r - — sinh r sin r 


2 4 

Cg = — cosh 2r cos 2r - 2 + - cosh r cos r 


C 3 " 


C 4- 


2 2 4 

— cos 2r - ■=• cosh 2r + ■=■ cosh r cos r 
5 5 5 


4 4 A. 

— cos 2r + — cosh 2r + ^ sinh r sin r 


(B36) 


2 12 
C5 = — sinh 2r sin 2r + — sinh r sin r 


2 12 
Cg = — cosh 2r cos 2r + 2 - cosh r cos r 


J 


a 

H 

11 

cosher 

cos 2 r 

sinh r 

d 2 = 

sihh 3 r 

sin 3 r 


d 3 = 

rz 

cosh°r 

cos 3 r 


D 4 = 

sinh^r 

sin^r 

cosh r 

D 5 “ 

cosh r 

cos^r 

sinh^ 

D6 = 

cosher 

cos r 

sin 2 r 

% = 

sinh^r 

sin r 

cos^r 

a 

CD 

11 

cosh r 

sin 3 r 

sinh r 

a 

to 

11 

sinh r 

cos r 


*10 = 

cosh r 

sin r 



(B37) 
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2 8 
E-j_ = ^ cosh 2r cos 2r + 2 - ^ cosh r cos r 

2 8 
Eg = j si 11 ! 1 2r sin 2r - — sinh r sin r 

2 2 8 

E„ = _ — cos 2r + — cosh 2r - ■=• sinh r sin r 

J 5 o o 

4 4 8 

E4. = - ■= cos 2r - •g- cosh 2r + — cosh r cos r 

2 8 

Eg = - -g sinh 2r sin 2r + — sinh r sin r 

2 8 
Eg = - g cosh 2r cos 2r + 2 - — cosh r cos r 

hj = cosh r cos r ^ 

hg = sinh r sin r 
h 3 = cosh r sinh r(l - 2 cos^) 
h 4 = - sin r cos r(2 cosher - l) 

L-j_ « (cosher + cos^r - l)^ 

Lg = (cosher - cos^r)^ (cosher + cos^r - l) 


V 


■> 

► 


M^_ = (cosher + cos^r - l)^ > 

Mg = (cosher - cos^r)^ 

M3 = (cosher + cos^r - l) (cosher - cos^r) > 

O O 

= cosh r sin r 

P p 

Mg ■«» sin^r cos^r 


(B38) 


(B39) 


(B40) 


(B41) 
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Ql = ^ cosh. 2r sin 2r 


^ cosh r sin r 


*\ 


2 

0,2 = ^ sinh 2r cos 2r 


sinh r cos r 


Q3 = ■=• sinh 2r cos 2r - 


■=• cosh r sin r 

D 


2 4 

Q4. = - ■£=■ cosh 2r sin 2r + — sinh r cos r 
o 5 


2 6 2 
Q5 = - g sin 2r + — sinh r cos r - — cosh r sin r 


06 = - — sin 2r + — cosh r sin r + — sinh r cos r 
5 5 5 

2 6 2 
Q7 = - g- sinh 2r + — cosh r sin r - — sinh r cob r 


> (B42) 


4 6 2 

Qq = -g sinh 2r - — sinh r cos r - — cosh r sin r 


J 


Nusselt Number Constants 

‘The constants appearing in equations (83) to (84) are given 
explicitly as 

O3. - (Ba-Bj.) sinh r cos r - (B 1 -tB 2 ) costl r sin r + 8(B 4 -B 3 )sinh 2r cos 2r 

- 8(B 3 +B 4 )cosh 2r sin 2r - 4B 5 (sinh 2r - sin 2r) (B43) 

Q 2 = (By -Bg) sinh r cos r - (Bg+B 7 )cosh r sin r + 8(B 3 +B A )sinh 2r cos 2r 

+ 8 (B^ -ll^) c °sh 2r sin 2r (B44) 

^3 = (B 9 -B 8 )sinh r cos r - (Bg-tBg)cosh r sin r + 8(B 11 -B 10 )sinh 2r cos 2r 

- 6(B 10 +B 11 )cosh 2r sin 2r + 4B 12 (sinh 2r + sin 2r) 


(B45) 
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ft 4 = (B^-B^jJsinh. r cos r - (B-jjg-tBp^cosh r sin r 
+ 8(B3-B^)sinh 2r cos 2r + S^g+B^cosh 2r sin 2r 

- 4Bg(slnh 2r - sin 2r) (B46) 

^ 5 = (B^g-B^gjcosh r cos r - (B^g+B-^sinh r sin r 

+ 8(B 18 -B 1 y)cosh 2r cos-2r - 8(B3_7+B2_g)sinh ^r s * n 

- 4B^g cos 2r + 4©20 cosh 2r (B47) 

^6 = ( B 22 " B 21^ cosi:1 r cos r “ ( B 21 + B 22^ sinl1 r sin r 

+ 8(B 1 74B^g)cosh 2r cos 2r - 8(B-)_7-B-i_g)Binh 2r sin 2r 

- 4B 8 o cos 2r - 43^9 cosh 2r (B48) 
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TABLE I- - ZEROTH OEEER UNIVERSAL FUNCTIONS 


M 


(a) Ra = 10 (b) Ra = 10 2 


n 

V 00 

V 01 

V 02 

v” 

01 

0 

0.8021 

0 

-0.3571 

0 

a 

.8835 

-.01293 

-.3536 

.07812 

.2 

.8876 

-.02504 

-.3431 

.1563 

.3 

.8946 

-.03570 

-.3256 

.2347 

.4 

.9041 

-.04387 

-.3011 

.3132 

.5 

.9159 

-.04902 

-.2694 

.3921 

.6 

.9298 

-.05028 

-.2304 

.4713 

.7 

.9457 

-.04676 

-.1842 

.5507 

.8 

.9628 

-.03771 

-.1304 

.6304 

.9 

.9811 

-.02237 

-.06901 

.7104 

1.0 

1.0 

0 

0 

.7906 


T] 

o 

o 

t* 

V 01 

Y 02 

v” 

T 01 

0 

0.36002 

0 

-0.59713 

0 

.1 

.36748 

-.03828 

-.59260 

.2220 

.2 

.38973 

-.07435 

-.57890 

.4463 

.3 

.42644 

-.10594 

- .55537 

.6753 

.4 

.47700 

-.13079 

-.52119 

.9108 

.5 

.54057 

-.14651 

- .4751 

1.1545 

.6 

.61597 

-.15069 

-.4154 

1.4073 

.7 

. 70173 

-.14079 

-.3402 

1.6694 

.8 

.79597 

-.11419 

-.2476 

1.9402 

.9 

.89632 

-.06818 

-.1350 

2.2180 

1.0 

1.000 

0 

0 

2.5000 


(o) Ra = 1600 


(dj Ea = 10 4 


1 

<1 

o 

o 

V 01 

y 02 



V 01 

r i 

v 00 

V 01 

v 02 



T 01 

0 

-0.1342 

0 

-0.1672 

0 

0 

-0.053806 

0 

0.022320 

0 

.1 

-.1258 

-.06668 

-.1735 

.09850 

.1 

-.056456 

-.012229 

.015544 

-.6621 

.2 

-.09985 

-.1322 

-.1929 

.2639 

.2 

-.062696 

-.031015 

-.0054315 

-1.2444 

.3 

-.05462 

-.1951 

-.2217 

.5616 

.3 

-.067254 

-.062081 

-.042058 

-1.6264 

.4 

-.01291 

-.2522 

-.2558 

1.053 

.4 

-.060949 

-.10908 

-.095256 

-1.6120 

.5 

.1061 

-.2986 

-.2886 

1.798 

.5 

-.030510 

-.17164 

-.16322 

-.9074 

.6 

.2279 

-.3266 

-.3105 

2 .840 

.6 

.040867 

-.24238 

-.23796 

.8741 

.7 

.3806 

-.3258 

-.3092 


.7 

.17153 


-.30121 

4.1660 

.8 

.5638 

- .2828 

-.2696 

5.897 

.8 

.37653 

-.32059 

-.32004 

9.3285 

.9 

.7745 

-.1807 

-.1736 

7.862 

.9 

.65964 

-.24313 

-.24306 

16.449 

1.0 

1.000 

0 

0 

10.00 

1.0 

1.0000 

0 

0 
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(b) IU - 10 2 



(d) H» - 10* 
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Figure 1. - Schematic sketch of configuration considered. 
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(a) m «* -1. 

Figure 2. - Dimensionless velocity distributions for various 
beat-source parameters with Ea = 10, = 10, C = -1. 



■V ■BPHHI 



y = Y/d 

(b) m = 1. 


Figure 2. - Continued. Dimensionless velocity distributions 
for various heat-source parameters with Ea « 10, K. => 10, 
C = -1. A 




Dimensionless velocity 


HACA TN 3141 


47 



y = Y/d 


(c) m = 2. 

Figure 2. - Concluded. Dimensionless velocity distributions 
for various heat-source parameters with Ra = 10, =» 10, 



Dimensionless velocity 
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(a) m - -1. 


Figure 3. - Dimensionless velocity distributions for various 
heat-source parameters with Ra = 10, K A - 10, C = -1. 
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Figure 3. - Continued. Dimensionless velocity distributions 
for various heat-source parameters with Ra = 10 , K. = 10. 
C = -1. 




Dimensionless velocity 
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y = T/a 


(a) m = -1. 

Figure 4. - Dimensionless velocity distributions for various 
beat-source parameters with Ra = 1600, = 10, C = -1. 



Figure 4. - Continued. Dimensionless velocity distributions 
for various heat-source parameters with Ra = 1600, K* =* 10 








Dimensionless velocity 
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co 

i 


pq 


o 



y = Y/a 

(a) m = -1. 


Figure 6 . - Dimensionless temperature distributions for vari- 
ous heat-source parameters with Ra = 10, - 10, C =■ -1. 
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7 = Y/d 

(b) is 1. 

Figure 6 . - Continued. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 10, 
k a = 10, C = -1. 
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(c) m =* 2. 

Figure 6. - Concluded. Dimensionless temperature distribu- 
tions for various beat-source parameters vith Ra « 10, 

K k - 10, C - -1. 
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y = Y/a 

(b) m b 1. 


Figure 7. - Continued. Dimensionless temperature distribu- 
tions for various beat-source parameters with Ra = 10, 
k a = 10, C = -1. 












Dimensionless temperature difference 


KACA 031 3141 


63 



7 = Y/d 
(a) m a -1 . 


Figure 8. - Dimensionless temperature distributions for vari- 
ous heat-source parameters with Ra = 1600, = 10, C * -1. 
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y = Y/a 

(c) m = 2. 

Figure 8. - Concluded. Dimensionless temperature distribu- 
tions for various heat- source parameters with Ra = 1600, 
k a = 10, C = -1. 


Dimensionless temperature difference 
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y = T/a 

(a) m = -1. 

Figure 9. - Dimensionless temperature distributions for vari- 
ous heat-source parameters with Ra = 10‘S = 10 > C ** _1 * 








HA.CA ra 3141 



Figure 10. - Dimensionless velocities for specisl case of 
M = 0, m = 1, «. 10, sad a = 10. (Note: For Ra - 10, 

the numerical solution u coincides with, given curve.} 
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Figure 11. - Dimensionless temperatures for special case of 
M = 0, m = 1, = 10, and a = 10. (Note: For Ra. =* 10, 

the numerical solution t coincides with given curve.) 
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log Ra 

Figure 13. - Nusselt numbers for wall 
at Y = d and K A » 10. 
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s 

o 

to 



0 .2 .4 .6 .8 1.0 


Figure 14. - Ratio of actual temperature difference to that 
for pure conduction for special case where M = 0 and m = 1. 






log Ra 


Figure 15. - Ratio of center-to-wall temperature difference 
to that for pure conduction for special care where M = 0 
and i*l. (Note: Dashed curve is plotted against yRa.) 
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